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1. Introduction 


Theta pinch systems occupy an important place in studies of dense 
high temperature plasma. Here, detailed theoretical study of tbota 
pinch, because of the necessity of tauing account of many nonlinear 
processes (gas dynamic motion, iragnotic field diffusion. Joule heating 

by temperature dependent conduction, electron and ion heat conduc- 
tivity, bulk energy losses and other effects), is impossible v\jthout the 
use of numerical calculations by computer. 

A two temperature nature is characteristic of theta pinches (and, 
generally, for fast, high current discharges). The ion component of 
the plasma in theta pinches is heated by rapid radial compression and 
subsequent slow final compression in a magnetic field which Increases 
over time. The electron temperature changes, both due to compression, 
and as a result of Joule heating. Since the relaxation time between 
the electrons and ions, as a rule, is greater than or comparable with 
the characteristic time of the process, the plasma components can have 
different temperatures: T^?^T^. 

Experiments shov; Xl~2] , that "end" losses (escape of mass, impulse 
and energy through the ends of the plasma pinch) also can significantly 
affect the dynamics of the process in theta pinches. 

To study a theta pinch with end losses, formulation of at least a /S_ 
two dimensional problem is necessary. However, fer a qualitative analy- 
sis of the processes, it also is advissible to consider a one dimensional 
problem, which models the end losses, with bulk escape of mass, impulse, 
and energy . 

^Numbers in the margin indicate pagination in the foreign text. 
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Sot.'<c' reaults of numerical calculations of theta pinch problems are 
presented in this study. The physical pi'occssos in theta pinch systems 
are considered in a cno d.imensJonal two temperature magnetohydrodynamic 
approximate on, Vfith allowance for end losses by longitudinal heat 
eenUuotivity . The numerical calculations are compared with oxper.1mcr.ts 
conducted earlier at SKTI £33 . The comparison gives satisfactory results 
for a number of parameters. However, there also are a number of dis- 
agreements of the results and experimental data. This m.akes further re- 
finement of the physical and mathematical formulation of the problem 
necessary. The results of preliminary numerical calculations of the 
two meter theta pinch built at SFTI are presented. 

2. Fo rmula tion o f the Problem 

Electrical capacitance C, with initial voltage Uq on it, is discharged 
in a coll, a circular metal cylinder cut along the generatrix. The 
cylinder length is and the radius, r^j. A discharge chamber is placed 
inside the cylinder. It is filled with a hot, conducting deuterium 
plasma v/ith Initial temperature Tg°=T^^. It is assumed that the chamber 
walls are electrically nonconducting and that the conductivity in the 
annular region between the plasma and the coil (vacuum) is sere. We 
will consider the process of compression and heating the plasma pinch with 
a magnetic field, on the assumption of axial symmetry, in a tv;o tem- 
perature ma:gnetohydrodynaraic approximation, with allov;ance for electron 
and ion heat conductivity. Joule heating and bulk energy losses. /J_ 

The corresponding system of inagnetchydrodynaml c equations, with 
allowance for the effects listea above, in Lagrangian mass coordinates, 
has the form; 
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Here, t la the tjme, r Is the Eu?.er cocrdinswte, m(dm«pzdr) Is the 
Lagrangian mass coordinate, the Lagrc-ngdan time derivative, v is the 
plasma velocit^y, p Is the donsitj. , ty, T„ end are the specific 
Internal onerp;y and temperature of the electrons and ions, respectively, 
and Pj are the electron and ion pressure, and are the radial 
heat flows due to electron and ion heat conductivity, k is the coeffi- 
cient of heat conductivity, is the o: noire r: -ion relaxation, Q,^ is 
the intrinsic plasma bremsstrahlung, and are the axial and 
azimuthal components of the magnetic field strength, respectively, a is 
the conductivity and c is the velocity of light. 


Functions W^. express the end losses of energy in the plasma pinch 
by longitudinal (directed along the pinch axis) component of the electron 
heat conductivity flow. It is expressed by the quantity 





=-x- 


<4- i. / 

/<f C t . if'' 


(9) 


whei’e K-K(Tg) is the temperature dependent Goex'ficlent of heat conductiv- 
ity and A is the length of the plasma pinch. 


We .shall consider the quantities v/hlch express the physical proper- 
ties of the deuterium plasma, or. the assumption of complete ionization, 
with the use of the formulas presented in ^,5 for their calculation: 
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In formulas (9)~(13), t-fo qua,ntli;y A Is tho CouXamb logarithm, the 
temperature is moasured in olectrcn volts, and the remaining quantitltes, 
in the CGf system. 


V/e will consider the equation of state of an ideal gas to bo valid: 

' (16) 


System of equations (l)-(l^!) is solved In the region 0< m< M(0<r<rQ), 
where r^ is the radius of the discharge chamber, M is the plasma masi 
in one radian, in. units of length of the plasma pinch. 
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(17) 


With m=0(r"0), the. synmetry conditions 
are assigned. With m=K(r-rQ), 


v,l:ire discharge current I(M,t) is determined from the electric circuit 
equations 
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with the initial conditions 




U(Az,/?y.i7^ 


( 21 ) 


Here, L, R and C are the inductance, resistance and capacitance in the 
external circuit, respectively, and is the Initial voltage. The 
problem formulated above was solved, by moans of difference methods 
developed at the Institute of Applied Mathematics, Academy of Sciences 
USSR, desc 2 -ibod in detail in [6-10] . We present some results of numerl 
cal solution of the problem formulated above. 
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3. Res u ltn oi‘ Num e rical Eyper3ments for “Small" Theta Plnoh 


A series of calculations v?as carried cut for the so called "sniall" 

theta pinch system, v-’hich. vms studied experimentally at SPTI 3 . The 

initial paranot'.rs vfore as follows: electrical capacitance G^lPOyP, 

initial voltage U~33kV, ohmic resistance of circuit R«1U“3 ohm, inductance 

L=7'10"^ At the initial time t=C, it is Etssumied that the deuterium 

plasma is heated to temperature T.'^'T. °=5eV, initial magnetic field 
0 © 1 0 — Q 

=0, chamber radius rn=3.5 cm, initial density p =2.34*10 ^g/enr 
(n =7'10' cm . The length of the plasma pinch along the z axis is 
jj=21cm. 


The changes of the outer boundaries of the piasma pinch r and dis- 
charge current I in the circuiiu vs. time during a half period are shown 
in Fig. 1. After the first maximum compression, which is reached at 
t=0,3 psec, some r’adial oscillations of the plasma are observed, v/hich 
die out completely at time t=1.4 psec. Current amplitude 1=1.16KA, and 
the discharge half period T/2=10.5 psec. 


The average \’’alues of temperatures T^ 
over the entire plasma pinch mass M, l.e.. 


and T and density p (averages 
.if J 7} e/rn etc. ) VS . 


time and number of neutrons 


C 


at time t . 


The maximum average temperatures are reached before the time of the dis- 
charge current: maximum (l'gjj.a_x~ 119eV), at t=4 psec and 440eV, 

at t=2.2 psec. This is connected primarily with the Icngitudlnal electron 
heat conductivity and, to a somewhat lesser extent, with the transverse 
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ion and electron heat ccnductlvltles . The maximum density p=2.9*10 ' 
g/cm” (p/p =124) is reached near the magnetic field maxjrrmm. The tem- 
perature and density fluctuations v;hich arise in the initial stage of the 
process. Just like the inertial oscillations of the plasma pinch (see 
Fig. 1), die out over timer. There ic a significant "breakaway" of the 
ion component of the temperature from the electron component, which /II 

reaches the highest values at the time of the maximum of T^^. The 
basic source of heating of the ion component is the v7ork of the compres- 
sion forces of the plasir.a, vrlth the iricreasing magnetic field. The ions 
exchange energy by collisions with electrons, part of which then is lest 


from the system by ler.p;j tudlna], electron heat conductivity (ste Fq. 3~^)* 
After the maicimum, the rate of increase in the maij;netic field is so 
small that the rate of escape of energy becomes pi'eater than the rate of 
supply of energy (sec Flh* < H'lorufore, toward the end of the half 
period, T> relaxes to the value of T . The number of neutrons per dis- 

•I- p* © 

chargej is the quortity K=0.6-10^ neutron/discharge. 


Profiles of mass Lagrnngian coordinate m as a function of temperature 
T^ (solid line), T^ (dashed line), density p (line v/ith triangles) and 
magnetic field strength H (dot-dashed line) at various moments of time, 
indicated by the corresponding figures on the oui^ves (time is presented 
in microseconds) are given in Pig. 3, In the initial stage of the 
process, the magnetic field forms a skin near the outer boundary of the 
plasma. The azimuthal current density and Joule heating rate are at 
a maximum in this region. Because of this, the electron temperature 
is higher than the ion temperature on the surface of the plasma pinch 
and lower than T^^ in the remaining part of the system. In time, the 
temperature equalizes quite rapidly over the mass of the pinch, and 
we have T^^>Tg everywhere. The radial distribution of the shock wave 
and its accumulation on the pinch axis are followed from tne density 
profiles. After the inertial oscillations of the plasma die out, the 
density distribution has a bell shape. 


The energy values which make up the balance of the total energy o' 
the system vs. time are presented in Fig. 5: 

C -t V •> Hi-d. -fJ2 "* S ~ ^ 
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elect roir.egnt 1 1 c energy in capacitance; 
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energy Irs inductance; 
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energy converted to Joule heat of external 
resistance (cii’’cult losses); 

^ magnetic field energy, kinetic energy and 

t 

^ total internal energy of system, respectively, 
at given bime; 






t 



enorfiy loDoes by lonsitudlnal olectron heat 
conductivity and due to plasma bremcotrahlungj 


3^ and 


initial circuit electromagnetic energy and 
plasma internal energy. 


It follows from the calculations that the efficiency of conversion 
of the electrical energy of the external circuit to internal energy, as 
a function of the initial conditions, is in the 0.0014-0.009 range. The 
energy losses due to longitudinal heat conductivity turn out to be 
comparable to the internal energy, while the energy losses due to brems- 
strahlung of a pure deuterium plasma do not play a significant part in 
the energy balance. The energy values which make up the Internal energy 
balance of the system vs. time are presented in Pig. 6: 

• (23) 


here, 
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0 0 


t- A/ 
r 




1/ am 




work of compression j 


/1 3 


o 'c 


Joule heat of system. 

It is evident that the work of compression plays a considerably greater 
part than Joule heat- 


Plg. 7 and Table 1 illustrate the effect of longitudinal electron 
(W^) and radial ion (W^) heat conductivity on the processes under con- 
sideration. The average ion temperature vs. time for various versions 
of the calculation is presented in Pig. 7, The case W_?^0, VI. /O is 
plotted as a solid line, W^=0, W^=0, as a dashed line, W^^SO, as 

a dot-dash line and W^=0, W^ao, as a dash-double tick line. The maximum 
plasma temperature and density, the number of neutrons per discharge 
and the energy conversion efficiency for the abovementloned versions of 
the calculations are presented in Table 1. It is evident that both 
current and current significantly affect the plasma parameters. 
However, the energy losses by longitudinal heat conductivity play a 
decisive role. Allowance for leads to a considerable reduction in 
the ion and electron temperatures. The efficiency of conversion of elec- 
tromagnetic energy to internal plasma energy is reduced to half, in 
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this case. HerOj the plaoma concentration chanseo comparatively little. 
The moot oonoltive parameter lo the neutron yield, which, depending on 
the calculation conditions, changes by three orders of magnitude. 

Currents and affect both the amplitude of the quantities and, also, 
these quantities vs. time (see Fig. 7). 


TABLE 1. 


Galculfltion 
Quan- ■ • version 
titles 

14'^' If 


! 1 

1 

• * i 

505 

7 j ? 

835 i 

i 


119 

II2 

1008 

903 j 

1 

.a' o' 

124 

III 

76 

75,5 ' 

I 

1 



3.4*10® 

2.5*10® i I.3*I0'^ 

; - - f 


0.0014 

0.00154 

0.0035 

0.0037 1 

1 

.lO-^secJ 

2.2 

2.42 

5.27 

1 

5.27 S 
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4. Comparison with Experiment 


Comparison of the calculations carried out for the "small" system, /I 6 
with allowance for the longitudinal electron and transverse ion heat 
conductivities, with the experimental results £3] shows that there is 
quantitative coincidence of the parameters (Table 2). This model des- 
cribes well the electrical engineering parameters of the discharge 
circuit (I, T) and the initial stage of the process. Including the rapid 
compression and subsequent inertial plasma oscillations. The calculated 
and measured electron temperatures are consistent. 

At the same time, significantly understated ion temperatures T^ and 
neutron yields and overstated plasma concentrations were obtained in the 
calculations. It does not correspond to experience at the stage of slow, 
adiabatic final compression of the plasma and time changes of the plasma 
pinch radius r, T^^ and N. 
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Xn thft experiment , r docreasoo monotonlcally , and and N reach 
maxlrimra valuoc Jn the maonefclc field maximum. In the calculations, r 
is practically constant ever time, but and N reach maximum values in the 
first quarter of the period. The cause of this discrepancy may be 
particle Icosoo through the ends of the system. The moot severe end 
losses should shew up in short pinches. Qualitative analysis shows 
that the escape of mass may result in additional heating of the particles 
remaining in the system, due to the adiabatic final compression of 
the plasma, 

TABLE 2. 


' I 

j - n° 

• v^] 

W dic- 
charge 
period 

[lO“^sec 

current 

i7 

[m] 

ist oscil- 
lation 
tirmo 

[l0“‘^sec 

1 

CO 

m 

^ e 
CeV] 

[oYj 

ynoutron 

discharge 

i Theory 


21 

1. 15 

0.3 

8.5*10^® 

119 

440 

0.6 ‘10® 

Experiment 



19+0.5 

1.3+0. 15 

0.3+0.05 

2*10% 

i.rxo^® 

100+20 

1200+200 

7' 10^2 *10® 


5. Results of Numerical Experiments for Two Meter Theta Pinch 


Similar calculations were carried out for the two meter theta pinch 
built at SPTI (ji=200 cm), experiments on which will start in the second 
half of 1979* The inital data were the following: C = 192.0yP, Ug = /17 

35kV, R = 3‘10"^ohm, L = 8.5’10"^H, T^° = T^° = 5eV. 

The average electron and ion temperatures and concentrations vs. 
time, as well as the total number of neutrons at time t, for initial 
density p°=2 . 3^ ‘ 10“" g/cm^ (n®=7’ lO^^cm"^) , are presented in Pig. 8. 

Version (1) corresponds to the case, when there is no escape of heat by 
longitudinal electron heat conductivity (W„=0). In version (2), the 
coefficient of longitudinal electron heat conductivity for a deuterium 
plasma is used (z=l). Allowance for longitudinal electron heat conduc- 
tivity reduces the value of T^ from 1900eV to 470eV. A comparison of 
the calculation results for the "small" (see Pig. 2) and "large" systems 
shows that a tenfold Increase in pinch length a, results in a 3-^ times 

increase of T„ (this result is consistent with the approximate analytical 
© 

relationship obtained from analysis of the energy balance of the 

o 


9 


olocti’on component of the placma) . The electron temperature reachoo the 
maximum value in the initial otage of the disohareo, and it oubooquontly 
changes little. 

A Dories of studies [11^ 12] in recoiit years has been devoted to 
study of a method of limitation of end losses, by means of physical 
plugs (closures) . Physical closures on the ends of the theta pinch 
should both exclude losses of particles from the plasma, and reduce 
(limit) energy losses caused by longitudinal electron heat conductivity, 
by the formation of a high z plasma in the area of the plugs 

, 1 + 0 , 09 .^-'' 

( ‘^// i 0 ,a.ss“^ value of z should be selected, so that 

radiation losses do not exceed losses by longitudinal electron heat 
conductivity. Sstimates show that, at the optimum z=30, coefficient /18 
K|| decreases approximately eightfold. Numerical calculations for z^ao 
(B'lg. 0(3)) give an app>eciable increase in electron temperature over 
the z=l case (this result should he considered approximate, since it v/as 
obtained within the scope of a one dimensional mathematical model with 
bulk discharges). 

The maximum ion temperature depends little on Ky . Since the re- 
laxation time is comparable to the half period T/2, a considerable 

gap is retained between T^ and T^, during the entire half period. 

The radial distributions of functions p, T„, T^ , and J, for 

0 2q o " 

time t=7.1 sec, with initial density p =2. 3^ * 10”^g/cm'^, are presented 
in Pigs. 9 and 10. With k|| =0 (Pig. 9), an electron temperature peak 
forms at the boundary of the plasma pinch (T layer [13] ) • Inclusion of 
the longitudinal electron heat conductivity corresponding to z=l results 
in resorption of the T layer (Pig. 10). 

The authors thank A. A. Samarskiy, R.G. Salukvadze and S.P. 

Kurdyumov for Interest in the work and useful discussions, as well as 
L.N. Busurina, for participation in programing, conduct of the cal- 
culations and processing of the results, 
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